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SUMMARY 


The solid-blockage interference for a doublet on the tunnel axis 
and the boundary interference for lifting wings in circular, rectangular 
and two-dimensional slotted tunnels have been calculated by substituting 
an equivalent homogeneous boundary for the physical boundary of discrete 
slots. In the case of small wings, the interference calculated with the 
assumption of homogeneity has been found to be consistent with that cal- 
culated for the discrete slots for as few as four slots in a circular 
tunnel. 

Furthermore, available experimental results for blockage interfer- 
ence are consistent with the results of the present analysis. As a con- 
sequence of the assumption of homogeneity it is possible to express the 
interference of multislotted tunnels as a function of a single parameter 
which combines the effects of two physical variables: the ratio of open 

to total slotted wall perimeter and the number of slots. A curve is 
presented which permits the rapid evaluation of this parameter and numer 
ical results for lift and blockage interference are plotted against the 
parameter . 


INTRODUCTION 


Several investigators have found that in a wind tunnel with bound- 
aries which are partly open and partly closed, the boundary interference 
on the lift of a wing can be reduced nearly to zero . References 1 to 5 
deal with the case of a doublet in the center of the wind tunnel, while 
references 6, ’J, and 8 consider the effects of wing span. There are 
several reasons why it would be desirable to have a wind tunnel with 
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zero-lift correction. In the first place, the necessity of applying the 
corrections would, be eliminated, although this, in itself, is not of 
large importance as long as the necessary corrections are known. More 
important is the fact that the pressure distributions cannot be corrected, 
and the distortion of both spanwise and chordwisc pressure distributions, 
as a result of boundary interference, places a limitation on the size of 
model which can be tested in a given wind tunnel. 

The primary interest in partly open or slotted wind tunnels, how- 
ever, is connected with the very different problem of wind-tunnel choking 
at high subsonic Mach numbers, which is a result of solid blockage inter- 
ference and which places a very severe limitation on the permissible 
model size for testing in closed wind tunnels. The blockage correction 
for a circular slotted tunnel (a wind tunnel in which the open part of 
the boun dar y is distributed around the periphery in the form of several 
longitudinal slots in an otherwise solid boundary) has been considered 
in reference 9 > which indicates theoretically that slot conf igurations 
exist for which the blockage correction is greatly reduced. The experi- 
mental results that are also included in reference 9 show that slotted 
tunnels can be used for aerodynamic testing in the transonic speed range . 
The conventional closed and open tunnels are both unsatisfactory in this 
speed range, the closed tunnel because of choking and the open tunnel 
because of the excessive power requirements and the large boundary inter- 
ference. Much effort has since been expended in the experimental devel- 
opment of transonic slotted tunnels, and several large tunnels of this 
type are presently in operation (refs. 10 and 11) or in construction. 

A knowledge of the lift-interference corrections for the slotted tunnels, 
as well as the blockage corrections, is thus of present interest. 

The lift interference in slotted tunnels has been considered in 
reference 12. This work has been extended in reference 13 to include 
wings of finite span, and numerical results are presented for several 
configurations of practical interest . 

The calculations required to determine blockage and lift correc- 
tions, for a particular slotted-tunnel configuration, by the methods of 
references 9 and 13 are very laborious especially for tunnels with a 
large number of slots . In order to obtain a general solution for this 
problem, Dr. A. Busemann of the Langley Laboratory has suggested that the 
problems of both lift and blockage interference be treated mathematically 
from the standpoint of a homogeneous boundary, with the slot effect uni- 
formly distributed over the surface of the boundary. It was reasoned 
that at some distance from the boundary, in the region of the model, the 
flow would no longer show the effects of the individual slots, particu- 
larly if the actual boundary contained a large number of slots . Further- 
more, the wind tunnels which are now in use contain several slots. It 
is felt that this is desirable in order to increase the uniformity of the 
interference. This paper will treat the problem in the suggested manner. 
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beginning with the development of a suitable boundary condition to 
represent mathematically a homogeneous boundary which has the same flow 
characteristics , at a point in the flow sufficiently removed from the 
boundary, as the actual physical boundary of alternate open and closed 
portions of the wall.- Numerical results will be presented for circular 
tunnels with slots uniformly distributed around the circumference, for 
rectangular tunnels with uniformly distributed slots in the top and 
bottom walls, for rectangular tunnels with the slot distribution deter- 
mined by a transformation from a uniformly slotted circular tunnel, and 
for a two-dimensional tunnel. 

The results of this paper are derived on the basis of an incompress- 
ible potential flow. The subsonic linearized compressible-flow theory 
shows that wind-tunnel lift-correction factors are not affected by Mach 
number; therefore, the lift-correction factors presented in this paper 
should apply directly to subsonic compressible flows, at least within 
the range of applicability of the linearized theory. The effect of com- 
pressibility upon the blockage interference is to increase the axial- 
interference velocities with increasing Mach number in proportion to the 

factor ~ ~ r/o ( ref ' 9) • 

(l - M 2 ) 3/2 


SYMBOIS 

The symbols A^ B n , C n , D n , C^, C 2n , Cj n , C W A k , A st , B st , 
C st> D st represent series coefficients. 

A area 

b semiwidth of rectangular wind tunnel 

c nondimens ional restriction constant 

c chord length of a two-dimensional airfoil 

Cj section lift coefficient 

C cross-sectional area of tunnel 

C L lift coefficient 

d slot spacing 

h semiheight of rectangular or two-dimensional tunnel 
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stagnation pressure 
restriction constant 
doublet strength 
number of slots 
static pressure 

open ratio of slotted wall (ratio of slot width to slot spacing) 

radius of a circular tunnel 
wing area 

thickness of slotted wall 

x-component of additional velocity due to presence of a model 
in wind tunnel, U - Uq 

free-stream velocity at upstream infinity 

x-component of velocity at any point 

y-component of additional velocity 

component of additional velocity normal to surface of wall 

z-component of additional velocity 
complex velocity in Z-plane 
distances in Cartesian coordinate system 
complex plane, x' + iy 1 
circulation around a wing 

Wo Q 

correction factor due to lift, — 

u o bC L 


£ complex plane, | + Itj 

Cartesian coordinates in £ -plane 
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9 angle in polar coordinates 

A height-width ratio for rectangular tunnel 

A 1 constant appearing in transformation from rectangle to circle 

u esc — r^ 

2 ° 

cp incompressible-flow velocity potential 

p radial distance in polar coordinates 

p density of air in wind tunnel 

o ratio of vortex span (effective wing span) to diameter of cir- 

cular tunnel or width of rectangular tunnel 

Subscripts : 

0 free-stream conditions at upstream infinity 

1 due to model in free air 

2 due to presence of tunnel boundaries 

C circular tunnel 

n in direction normal to wall surface 

R rectangular tunnel 

n,k,s,t summation indices 
r due to a row of vortices 

s signifies additional term required to satisfy boundary condi- 
tion at slotted walls 

T total 


BOUNDARY CONDITIONS 


Wall with discrete slots .- Longitudinal and transverse cross- 
sectional views of a slotted wind tunnel are shown in figure 1. In the 
longitudinal view exaggerated streamlines are drawn to indicate how the 
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air in the vicinity of the wall flows out and in through the slots as it 
passes the model. The chamber surrounding the tunnel is maintained at a 
pressure equal to the free-stream static pressxire of the flow inside the 
tunnel. The difference between this pressure and the average local pres- 
sure of the flow just inside the wall causes the air to move out and in 
through the slots. The boundary conditions at the slotted wall for this 
flow configuration will now be considered. 

Consider a set of axes in Cartesian coordinates which are fixed with 
respect to a model in the wind tunnel, while the air flows by with veloc- 
ity components U, v, w. Furthermore, let U = u Q + u, where u Q is 

the free-stream velocity at infinity. The pressure at a point inside 
the tunnel at the wall is given by 

p = H - ^-( u 0 ^ + 2u Q u + + ^ + 

The free-stream pressure p Q is given by 


Po 



Assume now that by some means the pressure just outside the wall is main- 
tained at p Q ; then, the pressure difference across the wall is 


&P = P - Pq 


n 

p v 


2u q u + 


u 


v^ + w^) 


Next, assume that the relationships between the model size and shape and 
the distance to the wall are such that u, v, w are all much smaller 
than Uq at the wall. The relation between the pressure difference and 

the axial velocity inside the tunnel then becomes 


Ap ~ -PUqU (l) 

This is a known result of the small-disturbance theory. Note that the 
small-disturbance assumption is required not in the field of flow near 
the model but only at the walls . 

It is now required to find an expression which relates the pressure 
difference across the wall to the flow through the slotted wall. This 
expression, combined with equation (l), will establish the relationship 
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between the axial velocity just inside the slotted wall and the flow 
through the slots. In order to find such an expression it will first be 
necessary to study the energy in the neighborhood of the slots . This 
energy is basic to the problem, for the essential mechanism of the slotted 
wall is that kinetic energy is stored in the air which flows out through 
the slots. This energy is later returned to the flow when the air flows 
back through the slots into the tunnel. Consider a thin slotted wall in 
a field of flow with a uniform velocity normal to the wall at infinity. 
(See fig. 2.) Because the flow pattern is the same for each slot it is 
permissible to study a single channel such as the one in which approxi- 
mate streamlines have been sketched in figure 2. The kinetic energy 
enclosed in a region of this flow bounded by a transverse plane at -x Q f , 

by the "walls" of the channel, and by the slot is given by 


Kinetic energy = 




9 



The region of integration may be considered to consist of a surface of 
unit depth normal to the plane of the page which includes the dashed line 
shown in figure 2. The component of velocity normal to the closed por- 
tion of the slotted wall must be zero. With regard to the open portion, the 
disturbance potential will be assumed to be zero at the slots, in con- 
formity with previous papers and with the classical practice in treating 

interference in open tunnels. Because = 0 at the channel walls, 

on 

and because cp = 0 at the slot, these two regions contribute nothing to 
the integral. With regard to the transverse plane at -x Q 1 it is clear 

that if this plane is sufficiently far away from the slot the potential 
will have a value which is essentially constant in this plane. 

Also, 



the quantity flow. Thus, 


Kinetic energy 


| P<P-x„ ,v n A 


"o 
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In order to complete the evaluation of the kinetic energy it is 
necessary to determine the value of the potential cp_ x * • Let the coor- 
dinate origin be taken in the plane of the slotted wall at the center 
of one of the solid sections and consider the flow to the left of this 
slotted wall . With the assistance of reference 1 4 , the potential of 
this flow is found to be 


ep = 



pqp IL r» 

c c 2 X ° 


S=J 


k x 

CLe b cot 



( 2 ) 


where k s = d is slot spacing, a is slot width, and r Q = 

The plane at -x Q ! is sufficiently far from the wall that the last term, 
which falls off exponentially as x* becomes more negative, may be dis- 
regarded. Thus, 



The insertion of this 
energy equation results in 

Kinetic energy 


V n (| log e CSC | r Q + x Q 

value of cp t in the 
~ x o 


^ h PI * lo S e C5C | r o + 


previous kinetic- 
x o ') v r? A 


In the absence of the slotted wall, the kinetic energy per unit area of 
the flow inside the region of integration would he 


1 

2 


P x o 



Consequently, the portion of the total energy which may be regarded as 
being due to the presence of the slotted wall is i piv n ^ per unit wall 

area, where 


1 = | lo §e csc | r o 


( 3 ) 


Note that the quantity l has the dimension of length. 
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Equivalent homogeneous vail .- Consider, now, the flow field that 
would result in the region between -x Q ' and the wall if the slotted 

wall were replaced by an imaginary homogeneous wall of zero thickness 
through which potential flow is possible and which is characterized by 

the existence in the local flow of an energy A pZv n 2 per unit wall area 

associated with a local normal velocity v n . The velocity of this flow 
field, at -x Q ' , will be essentially the same as that of the slotted- 

wall flow field; furthermore, the total kinetic energies between the 
planes at -x Q ' and at the walls will be the same . For the purpose of 

calculating the flow to the left of -x Q ' , therefore, the homogeneous 
wall is equivalent to the slotted wall. 


The result which has been obtained by studying the relatively simple 
flow from a source at is that there is associated with the flow 

through the slots a kinetic energy that is a function of the dimensions 
of the slotted wall and of the average velocity normal to the wall, and 
that this energy may be considered to be concentrated at the plane of the 
wall for the purpose of determining the potential of the flow at a point 
sufficiently far removed from the wall. If the singularity is located 
at a finite distance from the wall, or if singularities of types other 
than sources are introduced, there will be a velocity component parallel 
to the wall in addition to the normal component v n . Application of the 

principle of superposition shows, however, that the energy which is asso- 
ciated with a given flow normal to the wall will not be affected by the 
presence of additional velocity components parallel to the surface of the 
wall. It is necessary, though, that the slot spacing be small enough so 
that the difference in flow through adjacent slots is small. Thus, the 
analysis will be applicable only to tunnels with several slots. With 
this qualification the homogeneous wall will be equivalent not only to 
the slotted wall of figure 2 but also to the slotted wind-tunnel wall of 
figure 1 (section A-A) insofar as its effect on the model is concerned. 


If the slotted wall of figure 1 is replaced by an equivalent homo- 


geneous wall, the energy per unit area at the wall is — pZv 2 and the 

2 ^ 


momentum associated With this energy is pZv n per unit wall area. An 

individual particle of air which follows the outer streamline shown in 
figure 1 flows out and in through the wall as it passes the model as a 
result of alternate outward and inward accelerations due to the pressure 
difference across the slotted wall. The direction of flow of the normal 
component of velocity does not, in general, correspond to the direction 
of the pressure difference across the wall. Instead, it is the direction 
of the normal acceleration which corresponds to the direction of the 
pressure difference. This pressure difference, which acts in a direction 
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normal to the wall surface must, in the potential flow, be equal to the 
rate of change of the momentum associated with the presence of the slots. 
Thus, 


Ap = 




Although the quantity 1 is constant in time at a given point on the 
wall, it is left under the differentiation sign in order not to exclude 
the possibility that the slot configuration may vary from point to point 
on the tunnel wall. The only restriction which is placed on the axial 
or transverse variations of l is that they must not be too rapid, 
because the equation for l has been derived on the basis of a two- 
dimensional flow and a uniformly slotted wall. The derivative is given 
by 

mM = &K) + u !;K) + v £K) * * £K) 

Since the flow is steady and the slot configuration is constant in time 
when referred to the fixed axes, J^(^ v n) = 0. If now ^-^Zv n ) and 

JL(jv n ) are assumed to be of the same order as (or of higher 

order) then the acceleration is given by 

£K) ” “o IrK) 


to the same order of approximation as was used in obtaining equation (l) . 
The pressure difference across the wall is thus related to the velocity 
through the wall in the following manner: 

= pu ° |^( Zv n) 


Equating (l) and (4) for the pressure difference results in 
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Let the potential of the flow be given by $ = -u^x + cp, where cp is 

the disturbance potential. In terms of this potential the preceding 
equation becomes 


dep 

Sx 


eix V Sn^ 


or 




Integrating this equation in the x-direction along the wall gives 

Am 

9+1 ^ equal to a constant . One of the boundary conditions which 

must be satisfied by the wind-tunnel flow Is that there be no disturbance 
at infinity upstream, so the constant must be zero. Thus, the boundary 
condition at the wall becomes 


cp + 


Sep 

Sn 


0 


( 5 ) 


At this point let it be emphasized that, although this condition must be 
satisfied everywhere on the boundary, all the symbols, including Z, 
refer to only local values of the quantities represented. Thus, Z has 
been left free to vary in both directions on the surface . The boun dar y 
condition may therefore be used in the study of wind tunnels in which 
the slot width is variable in the axial direction. Furthermore, adja- 
cent slots may be of different width or different spacing. 

Although only a plane wall was considered in determining the rela- 
tionship between the wall dimensions and the wall-restriction con- 
stant, Z, it is possible to show that the same relationship holds for 

a circular wall. Consider the transformation £ = e^ applied to the 
region between y' = d/2 and y' = -d/2 in the flow field shown in fig- 
ure 2 . ^ By using polar coordinates (r,0) in the f; -plane, there is obtained 

r = e and 9 = y 1 . Thus, the line x' = 0 transforms into an arc of a 
circle of radius 1. (See fig. 3 *) Because of the linear transformation 
between y 1 and 0, the open ratio of the wall is unchanged by the trans-. 
formation. The source of the flow in the Z-plane at x' = -» transforms 
in the plane into a source at r = 0. The line x 1 = oo transforms into 
the arc r = «>. The expression for the potential inside the wall is 
obtained by applying the transformation to equation 2. 


<P = v. 


n 


log. 


;c — r. 


- log,, r + 


s=l 


kc, 

C s e 3 


log. 


cos 


2jts0 


(r < 1) 


(6) 
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The potential of the source at the origin in the absence of the wall 
would be 9 = -v log e r. Thus, it follows that the additional poten- 
tial at a point inside the wall, due to the presence of the wall, is 
nearly v n ^ log e esc r Q . Consequently, the restriction constant l 
for a circular tunnel is obtained in the same mariner as for a tunnel 
with plane walls. The function log e esc r D is plotted in figure k. 


Consideration of equation ( 6 ) will give some idea of the degree of 
approximation involved in assuming the additional, effective energy due 
to the presence of the slots to be concentrated in a plane at the wall. 

As an example, calculations have been made for a circular tunnel with 

8 slots, for which p = 1, where p = esc | r Q . (This gives an open 

ratio of slightly over 0 . 09 , which is in the range of practical interest.) 
Figure 5 shows the potential given by equation ( 6 ) plotted against the 
radial distance from the center of the tunnel. The calculations were 

made along radial lines through the center line of a slot ^9 = and 

through the center line of a panel (9 = 0) . Also presented is a curve 
calculated from the approximate or homogeneous wall equation 


cp = ~ log e p - log e r (7) 

Inspection of figure 5 shows that the approximation is quite satisfactory 
for the flow in the central part of the tunnel (say r <0.6). For tun- 
nels with more than 8 slots, the region of validity of the approximation 
will be even larger . 


ANALYSIS 

Circular Wind Tunnel 


Lift interference .- Consider a lifting wing symmetrically located 
in a circular wind tunnel which has a homogeneous boundary through which 
potential flow is possible . Let the wing be represented by a single 
horseshoe vortex. The downwash at the wing will be determined, in the 
classical manner, by finding the downwash in a tunnel cross section far 
downstream due to a pair of vortices of opposite sign with circula- 
tion r/2, where f is the circulation of the horseshoe vortex in the 
tunnel. The disturbance potential in the plane will be taken as cp^ + 
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where cp^ is the potential of the vortices in a free field and qp^ is 

the interference potential due to the walls. All length dimensions will 
be made nondimens ional by dividing by the tunnel radius R. The boundary 
condition at the wall (p = 1 ) is then 

*1 + "2 + c (|r + tr) = 0 (8) 

where c = the nondimens ional restriction constant. 

R 

The interference potential cp ^ must satisfy the equation V^cpg = 0 

throughout the interior of the tunnel. In polar coordinates this equa- 
tion becomes 


d 2f P2 + 1 ^2 + _1_ ^2 = Q 
dp 2 p p 2 de 2 

By the method of separation of variables, the following family of solu- 
tions can be obtained: 


oo 


sin n 0 


^2n cos 


) ( c 3nP n + C 4n p- n ) 


n =0 


Since cp 2 must be finite everywhere inside the tunnel, it must be f in ite 
at p = 0 and, therefore, C^ n = 0. Let 0=0 be in the plane of the 

wing. Then, because of the symmetrical location of the wing, it is 
apparent that cp^e) = q> 2 (rt - 0 ) • The cosine terms, which do not satisfy 

this requirement, are omitted by setting = 0 (except when n = 0 ) . 

The solution may now be put in the form 

A n p n sin ne'j ( 9 ) 



The potential of the vortex pair in a free field is, in polar coor- 
dinates. 
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cp = - JL tan" 1 2pP S - in J? - 

1 ^ 1 p 2 - ; 


( 10a) 


where a is the nondimensional semispan of the vortex pair. Also, 


^Pq _ r 2a sin o(p^ + a^O 

ap 4 " (p 2 - o 2 ) 2 + ltp 2 o 2 sin 2 e 


(10b) 


In order to find the interference, equations (lOa) and (lOb) can be 
expanded in Fourier series at the position of the wall. These series, 
together with equation (9), can then he inserted in the boundary condi- 
tion, equation (8), in order to determine the constants in the dis- 


turbance potential. The expansions for qp^ and 
the form 


d<Pi 

dp 


at p = 1 are of 


uu 

cp, = -tan -1 — sln - = \ B sin 
r 1 1 - a 2 1- 


n0 


n=l 


( 11 ) 


4* 5U 
r 


2a 


(l + a 2 ) 


sin 0 




sin n0 


(l - a*0 + 4p^sin^0 n=l 


( 12 ) 


No constant terms appear in 'these expansions so A Q = 0. By substi- 
tuting (9), (ll), and ( 12) in (8) the following equation is obtained for 
each value of n: 


sin n0 + B n sin n0 + cnA n sin n0 + cC n sin n0 = 0 


Thus, 


An - " 


B 


n 


+ cC n 
+ cn 


( 13 ) 


After A n is found, there remains the problem of determining the 

vertical interference velocity, and from it the interference factor 5. 
The vertical interference velocity along the line G = 0 is given by 
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(14a) 


n=l 


The average interference velocity between the origin and the point (a, 0) 
is 


Pa 


Wo = - - 




4 rt 

0 n=1 


dp 


oo 

“2 * - V "' 1 


(14b) 


n=l 


The interference factor 8 as given in reference 15 can be written as 


8 = 


W2C 

UoSC L 


where 

C cross-sectional area of tunnel (rt for the tunnel of unit 

radius ) 

S wing area 

Cj^ lift coefficient 

u c tunnel free-stream velocity at upstream infinity 


The circulation is related to the lift by 


SC T = its£ 


Using this relation there is obtained 8 = -2— w 0 

4ar ^ 

from equation (l4b) 


Substituting for W 2 
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5 = 1 y E P + cC n a n -2 ( 15 ) 

l 6 Z — 1 + cn 

n=l 


Instead of actually making the expansions indicated in equations (ll) 
and (12), it is possible to infer the values of the constants B n and C n 

from the known corrections for open and closed circular wind tunnels. 

For a closed circular tunnel Silverstein and White (ref. 15) give the 
equation 


6 = 


l6a 2 


log. 


1 + cx*~ 

'6 o 

1 - ex 2 


Making use of the series expansion (ref. l6) 


log. 


x + 1 

X - 1 


00 

-r 

n=l 


2 

nx 11 


(n — 1,3,5, • • •) 


p *i 

the following equation may be written after setting a = 


00 

6 = 2 ^ a 11 " 2 (n = 1,3,5, • • •) 


n=l 


(16) 


But equation (15) gives, for c = 00 } 


6 = 


16 


00 

r c v 


n=l 


(17) 


Because equations (l6) and ( 17 ) must agree , and since they can only agree 
if the series coefficients are identical, term by term, then 


c n = 2 o n (n = 1 , 3 , 5 ,...) 

C n = 0 (n = 2,4,6, . . .) 


(18) 
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The correction for an open circular tunnel is simply the negative of 
equation (l6) . For this case (c = 0 ) equation (15) gives 


CO 



n=l 


Comparison of equations (16) and (19) shows that 


( 19 ) 


Bn = - 


2g n 

n 


^ = 0 


(n = 1 , 3 , 5 , • • •) 

(n = 2, A, 6, . . .) 


( 20 ) 


Substituting equations (l8) and (20) in (15) results in, finally. 



( n — 1 , 5 , 5 , •••) (21) 


IBis correction is plotted as a function of the nondimens ional restric- 
tion constant in figure 6 for wings of small span. 

For the circular tunnel the nondimens ional restriction constant is 

given by c = ^ l°g e M-, where d is the angle between two successive 

slot center lines . Since the number of slots N around the circumfer- 
ence of the tunnel is given by 2jt/d, then 


c = l°g e M- 
N c 


( 22 ) 


for a circular wind, tunnel. 

Solid-blockage interference .- Consider next the problem of solid- 
blockage interference in the circular wind tunnel. The solid body in 
the tunnel is represented by a doublet with axis alined with the axis of 
the tunnel cylinder, and located at the origin of coordinates on the axis 
of the tunnel. The flow potential is again represented by cp 1 + q> 2 , 

where cp 1 is the free-field potential of the doublet and cp 2 is the 

interference potential due to the presence of the homogeneous wall. The 
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boundary condition which must be satisfied at the wall (p — l) , for all 
values of 0 and all values of x, is again given by equation (8) . 
Laplace's equation in cylindrical coordinates, which must be satisfied 
by cp 2 , is 


d 2 cp 2 i &P2 _ 

dp 2 P^P dx 2 


in the case of circular symmetry about the longitudinal axis x. The 
solution of this equation which will be used for the problem under con- 
side rat ion is 



where I and I-|_ are the modified Bessel functions. 

The free-field potential of the doublet is cp.. = - r-rr- 

(x 2 + p 2 ) 

The value at the wall of this potential and its derivative is now expanded 
in Fourier series form. 
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3x 


(x^ + l ) 5/2 


-n 


k=0 


sin 


totx 

L 


(26) 


T 1 *■ s l tCf) (27) 

The functions Q q and. are defined by equations ( 25 ) and ( 27 ) . Upon 

substituting the preceding equations into equation ( 8 ), the boun dar y 
condition becomes 



3x 


( X 2 + J5/S 


sin 



For each A ^ there results an equation 


A k = 



Substituting this in the equation for the interference potential gives 
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If, now, the fundamental wave length of this expression L is allowed 
to approach infinity, the summation can be replaced by an integral. The 
necessary relationship can be obtained as follows: 


Since 

Thus, 



(See, also, ref. 17.) By using this relationship and setting q = — , 
equation (28) becomes 


92 = 


^ 00 

2 

Qo(q.) - cQ x (q) 

n 

Jo 

I Q (q) + cql 1 (q)_ 


i Q (qp) sln q* d( i 


(29) 


The remaining task is to evaluate Qq and If equations (25) 
and (27) are integrated by parts the results can be brought to the form 
(see ref. 18) 


Q 0 (q) = qK 0 (q) 
Q 1 (q) = q 2 K x (q) 
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The final result for <p is, then. 


<Po = 


poo 

qK 0 (q) - cq 2 K x (q) 

Jo 

IqU) + cql 1 (q) 


I 0 (qp)sin qx dq 


The interference velocity is 


U2 = " 


&P2 2 

^OO 

a 2 

cqK x (q) - K 0 (q) 

-?J 

q 

0 

cqi 1 ( q) + i G (q) 


qx dq 


(30) 


This equation is for a doublet of strength = 1 in a tunnel of 
unit radius . For the general case, at the origin. 


O oo 


U2 


m 


2« 2 r5 


JO 


cqK-^q) - K Q (q) 2 

~~i — \ — r q dq 

cql x (q) + I Q (q) 


(51) 


where R is the tunnel radius. The axial interference velocity is 


plotted as a function of 



in figure 7* 


This parameter is 


chosen because it results in an approximately linear variation of the 
blockage interference . The values were obtained by mechanical integra- 
tion of equation (31) • Fortunately, the value of the integrand converges 
rapidly toward zero with increasing q. At q = 5, for instance, the 
value of the integrand was l/2 percent or less of its maximum value in 
the calculations which have been made. 


Rectangular Wind Tunnel With Top and Bottom Walls Slotted 

Lift interference .- Consider a rectangular wind tunnel of semiwidth 
unity and semiheight A. Inside the tunnel is a vortex pair of semi- 
span a located with the span parallel to the width (y) axis of the 
tunnel and with the center on the center line of the tunnel. The verti- 
cal walls at the sides of the tunnel are closed, but the horizontal walls 
at the top and bottom, with nondimens ional restriction constant Z, are 
partially open. The boundary condition at the closed side walls (y = ±1) 
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is, of course, ^ = 0. This condition is satisfied by a horizontal row 

oy 

of vortices along the plane z = 0 (the reflected images, out to y = ±oo, 
added to the vortex pair inside the tunnel) . The potential of such a 
row is given by 


cp r = 


J -tan~^ 

4 T 


tanh — z 
2 


tan ®.(y - a) 


+ tan' 


-1 


tanh — z 
2 


tan ^ (y + a) 


(32) 


and the vertical derivative by 


_ r 

8z 8 


>in it (y - a) 


sin it(y + a) 


cosh itz - cos it(y - a) cosh itz - cos jt(y + a) 


(33) 


To the potential cp r another potential cp s is added in order to 

satisfy the boundary condition at the horizontal walls. The total poten- 
tial cp = qp + cp o must satisfy the conditions 
i r s 


At y = ±1 


for _ o 
dy 


(34) 


At z = A 


T t + 


l 

bz 


0 


(35) 


At z = -A 


9 t 



(36) 


Since qp r already satisfies equation (34), <p E must also satisfy (34). 
A solution of Laplace ' s equation which meets this requirement is 




UU 

cos mty^A n sinh njtz + B n cosh rutzj + 


A q z 


n=l 


(37) 
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OQ) 

Because cp r is an odd function of z and is an even function, it 

dz 

is necessary that <p s be an odd function of z in order that both equa- 
tions (35) and (36) may be satisfied. Consequently, B n = 0 , and 


p 

<P S = 2_ cos sin ^ 1012 + AqZ 


ckp, 


s __ r 


Sz 


rutA^ cos ruty cosh rutz + A^j 


\n=l 


(38) 


Now, let equations (32) and (33) be expanded in Fourier cosine series 
at the boundary z = A. Then, 


cp r = ^ Y2 ° n C ° S mty 

n =0 


(39a) 


8cp r 

Sz 


DO 

= ilZ 011 


cos n jcy 


(39b) 


n=l 


(it is apparent from the form of equation ( 33 ) that there can be no con- 
stant term in equation (39b).) If equations ( 38 ) and ( 39 ) are substi- 
tuted in equation ( 35 ) the following coefficient relationships are found: 


Thus, 


+ ^ 0 A + ^o “ ^ 

C n + 2D n + A^ sinh rut A + ZrutA n cosh njtA = 0 


(n = 0 ) 
(n 4 0) 


Ao = 
An = 



A + 2 

C n + ^n 

sinh rat A + 2 rat cosh rat A 


r (4o) 
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Let cp T = cp-^ + cp 2 where again cp x is the free-field potential of the 
vortex pair inside the tunnel and cp 2 is the interference potential. 
From the free-field potential of the vortex pair 


4n 


k-n _ y_-_£ 


y + a 


^ f Sz (y - a) 2 + z 2 (y + cr) 2 + z 2 


Thus, the interference velocity at any point (y,z) i£ 


4n _ « 
r ^ 2 


sin *(y-o) 


cosh itz - cos tf(y - a) cosh jtz 


sin rt(y + a) 

jtz - cos jt(y + a) 


E 

n=l 


y - a y + a 

cos mty cosh natz - A Q - — + 


(y - a ) r> + z 2 (y + a) 2 + z 2 

(U) 


The average interference velocity in the plane of the wing (z = 0) is 


W 2 . - I w 2 dy 


In order to evaluate this integral, obtain first 


ki t 

r 


v 2 ay - logg i - cos dz - Si - 2 y~ 
^ e 1 - cos it (y + a) / 


-y 


v a sin rmy 

n=l 


2A 0 y + log e 


(y + g ) 2 
(y - a ) 2 


As y approaches a the first and last terms approach minus and plus 
infinity, respectively. In order to eliminate the indeterminateness, 
the two terms are combined and the cosine is expanded in series form. 
Finally, the limit is taken as y approaches o. 
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lim iZ + 111 1 - cos A (y - q) 

y — >-a 1 - cos *(y + °) (y _ a ) 2 


lim + 0)2 

1 

H 

"S' 

1 

1 

C\J 

b 

1 

J>> 

CVJ 

K 

4 

r“ 

1 

If 

2 41 

y j .t 1 - COS it (y 4- a) 

_ (y - o) 2 _ 

W 2 ,2 


1 - cos 2jta 2 



By using this limit, the average interference velocity becomes 


w 2 



OO 


n-l 


sin mto + 


_C 2 o_ 
A + l 


( 42 ) 


For this rectangular tunnel 



so 


k-Tia 2 


Correction factors are plotted as a function of the restriction 
constant for several values of a in figure 8 (a) for a square tunnel 
(A = 1 ) and in figure 8 (b) for a tunnel with A = O. 5 . The nondimen- 

sional restriction constant is defined as c = where h is the semi- 

h 

height of the tunnel. For the tunnel under consideration then, c = — . 

A 

The value of c for a rectangular tunnel of height-width ratno A with 
N slots in each horizontal wall is given by 


00 

lQ ge( : — ) ~ ^ ^ 

\sin ttct/ L 

n=l 


sin rut a + 


(a + l) 


(43) 
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Solid-blockage interference .- Consider next the problem of solid- 
blockage interference in the rectangular wind tunnel with top and bottom 
walls slotted. The solid body in the tunnel is again represented by a 
doublet with axis alined with axis of the tunnel and located at the ori- 
gin of coordinates on the tunnel axis . The boundary condition at the 


closed side wall (y = ±1) is 



This condition is satisfied by 


a horizontal row of doublets with axes alined with the axis of the tun- 
nel and placed along the line x = 0, z = 0 at y = 0, y=±2, 

y = ±4, . . . y = + 00 . The potential of such a row is given by 


OO 



k=-oo 



X 

+ 2k) 2 + z 2 


+ x§ 5/2 


and the vertical derivative by 


(^ 5 ) 


c*Pr 

dz 



(46) 


The total potential cp T = cp r + cp s must satisfy equations (34), (35), 
and (36). As before, cp s must satisfy equation (34). A solution of 
Laplace's equation which meets this requirement is 



(47) 


As in the case of lift interference cpj. is an odd function of z 


and 


dz 


is an even function; therefore, cp s must be an odd function of z 


in order that both equations (55) a-ncL (5&) satisfied. Conse- 

quently, B s ^ =0, and 
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s=0 t=0 


A s ^. cosh 


f“ )2+ (rf 


cos srty sin 


(48) 


8cp 


s _ m 


6z kit 


3=0 t=0 


A st \ 

(s *) 2 4 - (- 

— ^ sinh 
L / 

z \ 

(sn ) 2 + (- 

tf 


\ 



\ 



cos sity sin x 
L 

(49) 


Now let equations (45) and (46) be expanded in Fourier cosine series 
at the boundary z = A. Then, 


00 00 


W Z=A = - ^ Y1 H c st cos s «y sin ~ 


s=0 t=0 


(50) 


’<*Pr 

.dz 


00 00 


z=A 


s=0 t=0 

where C s -j. and. D s ^ ane given by 

U pL 


]—■ ^ ^ D st cos s «y sin — 


r _ % 

^st “ T” 


'0 J 


I 


0 


x sin cos sjty 

Ju 

[Ty + 2k) 2 + A 2 + xf] 5/2 


dx dy 


(51) 


where 


1/2 

P = / 


(s = 0) 

(s 4 0 ) 


c 3t - k Mf) 


^0 


pi 


PL 


cos ^2^ cos sjty dx dy 


0 


f(y + 2k) 2 + A 2 + x 2 


(52) 
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D 


st 


4p 


m pl 


3 Ax sin cos s«y 


k=-co 


0 u 0 


(y + 2k) 2 + A 2 + x 2 ^ 


dx dy 


P 1 


_ |fP/rtWtn 

« v l A l 


u 


pL 


0 u 


oo 

Z p 

k=-» ^ 


A cos cos sny 


(y + 2 k ) 2 + A 2 + x 2 ] ^ 


dx dy 


(53) 


By substituting equations (48), (49), (50) , and (51) into equation (35), 


A s .(. can be written as 


^st ” ^st 


v st 


cosh 

A \ 

f s 2 /tj(\ 2 

(sjt) + hrj 

+ l J(srt) 2 + (- 

ffsinh 






1 




(5*0 


By substituting the equations ( 52 ), (55), and (54) into equation (48), 
the result is 





. f 1 

ta V" 1 

L | j 

‘-'0 w 

j IX COB 

^ cos Bird cot 

; cos sna 

L 

! [(,,♦ 2k) a 

0 L 

♦ ^ + e g 5/2 g a * ■*)* * >? ♦ 

cosh jx ^ 

|t-» a * (rf 

+ l \ (b*)^ + sinh 

H u, ' a * (t ) 2 


w_° 


(55) 


As L approaches 
becomes 


00 , the integral with respect to L in equation ( 55 ) 
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(V7) 
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LcA,Efe* ao g _+ * 2 ] . Ko r j (a t , k): > + A ajl 

^(a * 2k) 2 + X 2 J 

oohjxjq 2 + U*)^ + cA>jq 2 + (B*) 2 sinh Jx ^ 2 + ( 


J :-ot;h z\Jq 2 + (sk ) 2 cos 


3 Sl 


q 2 jqcA 2 Hi(X,s,q) - H (X,s,q)l 1-5 : 75 . 

L- 1 =T cosh z\fa + (s«)^COB qx coo s*y dq 

C(X,o,q) + cS(X,s,q) ^ 


K^a + 2k) 2 + A 2 ] 


COS STiOC da 


J(a + 2 k ) 2 + A 2 


’1 . 00 


i>y (a + 2 k)^ + A^J cos sjta da 


'0 k= 


C(A,s^q) = cosh A\jq 2 + (sit ) 2 


S(A,s,q) = A\jq 2 + (sit) 2 sinh A\Jq 2 + (sit) 2 

and where 2 has been replaced by cA . The interference velocity \x^ 

of a doublet symmetrically located in a tunnel with its axis alined with 
the axis of the tunnel is given by 

lip = -b u s - u-^ 

= . . A(„ . 

dx Sx\ r 1' 


dx + 4 * 


OO 

— y 

k* z — 


(y + 2k) 2 + z 2 - 2x 2 


|2 + z 2 + x 2]5/2 


(k 4 0) 
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This is for a doublet of strength m/4j( in a tunnel of unit semi- 
width. For the general case, at the origin. 


"2 = - ir 1 


v x /x=0 
y=o 
z=0 


m _1_ 

32«b3 Zy k 3 


(59) 


Equation (59) was solved numerically for the interference velocity U 2 
of a doublet in a square tunnel (A = 1) , and the ratio is plotted 

{ 


against 


c + 1 


in figure 9* 


Lift Interference in Rectangular Tunnel 
With All Sides Slotted 

Consider the lift interference for a small lifting wing symmetrically 
located in a rectangular wind tunnel of semiheight A and unit semiwidth 
and with all four sides slotted. The widths of the individual slots are 
determined by calculating the points on the walls of the rectangular tun- 
nel which correspond to the slot edges of a uniformly slotted circular 
tunnel, with the assistance of a transformation which maps the perimeter 
of a rectangle of semiheight A and unit semiwidth into the perimeter 
of a circle of unit radius. In the analysis of this particular problem 
the transverse cross-sectional plane of the rectangular tunnel, which has 
been previously the y,z-plane, will be taken as the complex Z-plane where 
Z = x + iy. The cross-sectional plane of the circular tunnel, which has 
previously been treated in terms of the polar coordinates (p,0) will here 
be taken as the complex £ -plane where £ = £ + irj . In reference 20 the 
transformation which maps a rectangle in the Z-plane into a unit circle 
in the £-plane is given as 


t = sn A'Z dn A'Z 
S cn A*Z 


(60) 


or 


p 2 = 1 - cn 2A*Z 
1 + cn 2A’Z 


(6l) 
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and A' is determined from 


K = K1= A' (62) 

2 2A 


where K and iK' are the quarter periods of the preceding Jacobian 
elliptic functions. A method for calculating A' and tables of the 
preceding Jacobian elliptic functions are also given in reference 20. 

The quadrantal perimeter of a square (A = l) and the quadrantal perimeter 
of a rectangle (A = l/2) were mapped into the quadrantal arc of a unit 
circle by using equation (6l) and the results are given in figure 10. 

The lift-correction factor for a small lifting wing symmetrically 
located in a rectangular wind tunnel may be written as 


Av, 


5 R ~ 


2R 


A (v 2R ) 


Z=0 


a R f 


a R f 


(63) 


since v^ ~ ( v 2r) r ° r a sma11 lifting wing. 

Z— 0 

The complex velocity in the Z-plane at a point on the x-axis is 


dW 

dZ 


-IVr = 


dW dC 
d£ dZ 


-xv c 


dZ 


or 


V R = V C 


dZ 


(64) 


The interference velocity for the rectangular tunnel is thus 


'2R 


= ( V 1C + v 2c)t| " V IR 


(65) 


In reference (21) d£/dZ is written as 


A 1 

dZ 


( 66 ) 
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at the origin. Evaluation of equation (lOb) at the origin gives 



(67) 


From the analysis of lift interference in a circular tunnel (eqs. (l4a), 
( 13 ), (l8), and (20)) 


V 2C = 


Fg C /c - 1 \ 
2* Ip + l) 


(68) 


at the origin. The vortex semispan in the circular tunnel a c is deter- 
mined from equation (6l) by setting Z = (whence, £ = a c j . The 

velocity v^ is given by 


v 


1R 


f 

2lfa R 


( 69 ) 


at the origin. This differs from v ic (eq. (67)) only because of the 

change in vortex span. If, now, equations ( 66 ), (67), ( 68 ), and (69) 
are substituted in equation ( 65 ) the result is 


V 2R = 



+ a r 


c + 1 J 



(70) 


Substituting this in equation ( 63 ) gives, for the lift-correction factor 
at the center of the vortex span. 


6r = 



c - 1 
c + 1 



(71) 


In the limit as 


0 this becomes 
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( 72 a) 


where m is the parameter (see ref- 20 ) of the elliptic functions of 
the transformation. The behavior of m is such that if A < 0.5 it i s 
sufficient to use 



When A = 1.0 (square tunnel), m = 0.5 and 


( 72 b) 


&R - 



c + 1 


= 0.274 




When A >2.0, the correction factor is approximately 


6 


R “ 



1 \ = n \(2 _ 1 ) 

c + 1/ l6\5 c + 1/ 


( 72 c) 


( 72 d) 


Equations ( 72 ) apply to wings of small span. 

The lift correction at the center of the wing has been calculated 
for a small wing in a square tunnel (A = l) and for a small wing in two 
rectangular tunnels (A = 0.5 and A = 2) ; the lift-correction factor 

is plotted against — — — in figure 11. 

c + 1 


Two-Dimensional Tunnel 

Lift interference .- Consider the problem of the lift interference 
in a two-dimensional wind tunnel with the top and bottom walls slotted. 
The wing in the tunnel will be represented by a two-dimensional lifting 
vortex located at the origin of the coordinate system, which is centered 
between the slotted walls . The potential of this vortex in a free field 
is 


©1 = — tan - -* - — 
T1 2* x 
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and the vertical derivative is 


^1 = r_ x 

2* x 2 + z 2 


The interference potential takes the form 


<p c 


00 

= ^ sin 222L sinh ^ 


A z 
o 


n=l 


( 73 ) 


There is no axial interference along the tunnel center line (z = 0) 
ckpo 

because — — = 0. The vertical interference at the origin is simply 


v 2 = 


&P2 

Bz 


= “Ao 


Inasmuch as the summation does not contribute to the wall-induced 
velocity at the origin, it is sufficient to study the term AqZ by 

itself . If the multiple-valued function for cp x is set eq ual to zero 
at the upstream infinity (z = -<») it is apparent that the Fourier series 
expansion of <Pq at the top wall will result in a constant term - — 

4 

plus an odd function of x. The expansion of ckpw^z will be an odd 
function with no constant term. Thus, if the constant terms in cp^, <p o, 
and S(p 2 /^ z are inserted in the boundary condition at the upper wall 
(z = h) 


<P. 


T 


+ 2^=0 

oz 


(74) 


- £ + V 1 + 2A o = ° 


A 


o 


X 1 

kh c + 1 


the result is 
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where 

by 


c = The vertical interference velocity at the origin is given 
h 


w 2 


= _ JL L_ 

kh c + 1 


Thus, the interference is found to be downwash at the origin which 
varies in magnitude from zero for a closed tunnel to a maximum for the 
open tunnel. Because 


r 


c Z u o c 

2 


the induced-flow angle at the origin is 


1 £ 

8 h c + 


where 

c = ^ lo g e csc | r o 

and 

d slot spacing on horizontal walls 

h half-height of tunnel 

r Q open ratio of the uniformly slotted horizontal walls 

Cj airfoil-section lift coefficient 

c airfoil chord 


(75) 


The induced angle is, of course, not constant along the chord of the 
airfoil; therefore, the wail-induced flow can be considered to have a 
certain curvature . This curvature has an effect similar to that of camber 
in an airfoil, and if the length of the chord of the airfoil is suffi- 
ciently great compared to the height of the tunnel, a curvature correction 
is required (ref. 3) . This curvature can be determined from a complete 
solution for cpg (eq. (73)) under the boundary condition of equation (7*0- 
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Solid-blockage interference . - Consider the problem of solid-blockage 
interference in a two-dimensional wind tunnel with top and bottom walls 
slotted. The solid body in the tunnel will be represented by a two- 
dimensional doublet with axis alined with the tunnel axis, and located 
at the origin of coordinates on the tunnel axis. The potential of this 
doublet in a free field is given by 


tpi = _ 


m 


^ x^ + 


(76) 


and the vertical derivative by 


dq>i 

dz 


m 


xz 


” (x a + z 2) ; 


(77) 


The interference potential must satisfy Laplace 1 s equation. Since 

qp^ is an odd function of x and an even function of z, qp£ must also 
be an odd function of x and an even function of z. Therefore, <P 2 
can be written as 


00 

cp 2 = ^ ^ A k cosh ^ sin 


krtx 


n=l 


(78) 


and the normal derivative by 


8cpg 

dz 


00 

= — Y A k sinh 
2 jt Z — L 


kjtz sin k«x 


n=l 


(79) 


The boundary conditions to be satisfied along the slotted walls (z = ±h) 


are 


at z = h 


<Pi + q >2 + 1 ^-(<p 1 + 92) = 0 


(80a) 


and 


at z 


-h 


q>i + q > 2 - 1 j^-(q>i + cp 2 ) = 0 


(80b) 
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These conditions are satisfied by expanding cp^ and into 

Fourier sine series in x along the slotted walls and then solving for 
the coefficients of cp 2 that satisfy these equations. Along the 

line z = h, equations (76) and (77) can be written as 


where 


and 


Vi - - £ 5! sln T 

k=l 


^>1 _ jlV" 

8z 2it Z — 


C v sin 


kjtx 


k=l 


PL 




uo 


(x^ + h^) 


sin dx 


PL 


hx 


0 U 2 -* 2 )’ 


sin dx 


A single integration by parts results in 


(81) 

(82) 


( 83 ) 


Ck = 


knh 


PL cos 

L_ 

x 2 + h 2 


dx 


( 84 ) 


By substituting equation (83) for and equation ( 84 -) for into 

equations (8l) and (82), respectively, and then solving equa- 
tions (80), (8l), and ( 82 ) for A k , the result is 


A k - 



cos 


kjtx 


- sin 


krtx 


dx 


X 2 + 


1 kfl 


sinh 


krth 

L 


(85) 
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By substituting equation (85) for A^ into equation (78), the result is 


^2 = “ o 



L /hZkjt krta . kna\ 

~ C ° S ~ - a Sin T~ ) 
V a 2 + h 2 / 

cosh felil + i iS3L sinh 


da 


cosh sin 

L L 


k=l 


( 86 ) 


Let q = £q. = and let L >■ ». Then equation (86) becomes 

Li Li 


n°° 


qp 2 = - ^ 


/hZkrt 

Jo v 


T cos qa - a sin qa \ 

Ida 

cosh qh + Zq sinh qh 


a 2 + h 2 


cosh qz sin qx dq 


60 


Zqe _< lL _ e -qh 


m 

2« | \cosh qh + Zq sinh qh, 

10 


cosh qz sin qx dq 


m 

2n 


e~^ h ( Zq - 1) 


cosh qh + Zq sinh qh 


cosh qz sin qx dq 


m 

2*h 


noo 


e“^-(cq - l) 


cosh q + cq sinh q 


UO ■- 


cosh sin dq 
h h 


( 87 ) 


where qh has been replaced by q and c — — . 

h 

ity 122 is 


The interference veloc- 
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^ ox 2r.h 2 


Ooo 


qe _c i(cq - l) 


cosh — cos — dq (88) 
cosh q + cq sinh q h h 


At the origin, equation (88) reduces to 


U2 = 


m 


2«h c 


noo 


qe~^(cq - l) 
cosh q + cq sinh q 


dq 


(89) 


0 


Equation ( 89 ) has been evaluated numerically, 
velocity function Ugh^/m is plotted against 


and the 




1 

c + 1 


axial- interference- 
in figure 12 . 


DISCUSSION 

Application of the Results of the Interference Analysis 


Consider the problem of determining, say, the correction factor 6 
for a particular wing in a square tunnel. The tunnel has four slots in 
each horizontal wall, with an open ratio of 0.164. The vertical walls 

are closed. From figure 4, log e esc | r Q = 1-37- After using equa- 
tion (44) to determine c, the value of — ^ — is computed to be 0.82. 

The wing has an effective span of one-half the tunnel width, and for this 
case, figure 8(a) gives the value of the correction factor as - 0 . 062 . 

If a semispan reflection-plane model of a wing of twice this span is 
tested, the effective height-width ratio of the tunnel will be 0-5 and 
the correction factor will be -0.14-6. 


Comparison With Previous Results 

A comparison of the results for the lift interference of circular 
tunnels with the results of reference 13 , in which the individual slots 
are considered, shows striking agreement for wings of small span (fig. 6) . 
In the case of a lifting doublet in the center of the tunnel the results 
agree within 2 percent for only 4 slots and even more closely for greater 
numbers of slots. This agreement would seem to be sufficient proof of 
the justification for the basic assumption underlying the present theory, 
that the wall can be considered to be homogeneous. This assumption is 
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not made in reference 13 . The theory of this reference indicates that 
for small numbers of slots, say 2 or k, the lift interference varies 
appreciably with the orientation of the slots. This variation, of course, 
cannot be predicted by the present theory. 

The results of the analysis of solid blockage interference are com- 
pared with the results of the analysis of reference 9 in figure 7 . The 
results of this reference are seen to straddle the curve obtained in the 
present analysis . The spread in the results of reference 9 may be due 
to the fact that the series used in the calculations converged very slowly 
in the case of ten slots . 

The present theory will now be compared with the experimental results 
of reference 9» For the 10-slot circular tunnel the present theory pre- 
dicts a small negative interference. The magnitude of this predicted 
interference is about l/l 3 that for a closed tunnel and about 1/3 that 
for an open tunnel. An inspection of figure 7 (ref. 9 ) shows that the 
measured interference in the slotted tunnel checks very closely with 
these values. This is a single-point comparison at the midpoint of the 
body, however, so it can not be regarded, in itself, as a complete veri- 
fication of the theory . In fact, this body was so large with respect to 
the tunnel diameter that the present theory, based on a doublet in the 
tunnel, could hardly be expected to apply. 

A much smaller body was tested in the 8 -slot octagonal tunnel of 
reference 9 • For this tunnel the present theory predicts the value of 
the axial- interference-velocity function to be -0.006. This is less than 
5 percent of the magnitude of the interference in a closed circular tun- 
nel. The predicted interference velocity ratio for the 1— - inch model 

3 

in the slotted tunnel is — • = -0.0002 for incompressible flow. Even 

when multiplied by the compressibility factor (1 - M 2)-3/ 2 , it is obvi- 
ous that this correction is within the usual experimental accuracy for 
Mach numbers within the range of application of the linearized theory. 

For instance, at M = 0.9j the axial interference velocity is approxi- 
mately 1/4 of 1 percent of the free-stream velocity. The experimental 
results in figure 12 of reference ^ show that the interference was negli- 
gible along the whole length of the body at all speeds up to a Mach num- 
ber of about 0.9, Thus, these experimental results are consistent with 
the present theory with regard to blockage interference . 


Application of the Concept of a Restriction Constant 

The concept of a restriction constant can be useful in other ways 
than in the direct solution of problems in slotted wall interference as 
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for instance, in providing a common basis of comparison for walls with 
different slot spacings. As an example of this type of problem, consider 
a circular tunnel with eight evenly spaced slots around the circumfer- 
ence, with an open ratio of 0.15- Suppose it is desired to construct 
another circular tunnel with the same interference characteristics as 
this tunnel but that it is also desired to increase the number of slots 

to 12. For the given 8 -slot tunnel c = | log e p = O. 365 . For the same 

value of c, the 12 -slot tunnel must have an open ratio of 0 . 071 - 

Using this concept, it is also possible to determine the effects of 
changing the shape of the slot cross section. With the realization that 
the theory depends upon an essentially potential flow in the slots, for 
instance, the question arises as to whether the sharp-edged slots could 
not be improved upon. It might be advisable to try smoother shapes in 
an attempt to obtain a closer approximation to potential flow in the 
slots . The restriction constant Z could be determined for any slot 
shape either analytically by calculation of the potential flow through a 
group of such slots in cascade or experimentally by a simple electrical 
analogy. Once Z has been determined the results of the present analysis 
could be used to predict the interference for such slots . 


Experimental Determination of the Restriction Constant 

If experimental data of sufficient accuracy to permit the evaluation 
of the lift- or blockage-correction factor for a particular slotted tun- 
nel are available, the effective value of — i — can be read directly 

c + 1 

from the curves presented in this paper. Once the value of this parameter 
is known, corrections may be computed for other wings in the same tunnel. 


Limitations of the Theory 

Some question arises as to the range of validity of the equation 
for the restriction constant of a slotted wall, inasmuch as the open ratio 
required to obtain a given value of c decreases rapidly as the number 
of slots increases. Since the normal mass flow remains constant with c, 
this requires a continuously increasing normal velocity in the slots 
which, although it is permissible in the assumed incompressible flow, is 
not possible in the actual physical flow. It appears, from these consid- 
erations, that the limit of applicability of the restriction-constant 
equation will be determined by the required normal flow velocities rather 
than directly by the smallness of the slots. This limitation should be 
considered in the application of the results of this analysis to the 
design of wind-tunnel test sections. 
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There is, also, some question regarding the use of potential-flow 
theory in the region near the slots . This question arises because the 
potential theory omits certain flow phenomena, in particular, the pres- 
ence of a mixing region at the boundary between the moving tunnel air 
and the essentially static air outside. The adequacy, for the purpose 
of interference computations, of the potential theory must be determined 
by experimental studies . 


The restriction constant has been determined by assuming the effec- 
tive free boundary (cp = 0) to be located along the line between the slot 
edges. If the boundary were actually located outside of this line, the 
slot restriction constant would be somewhat larger. If the slots had a 
finite depth t, instead of being essentially sharp-edged orifices, and 
if the effective free boundary were at the outside of the slots, the 

restriction constant would be approximately Z = — log u + Actu- 

Jt e r o 

all y> it is likely that in the case of slots of considerable depth, the 
location of the effective free boundary would be a function of the local 
outflow or inflow. Not only would this make any calculation of the inter- 
ference questionable but it would also present the distinct possibility 
that the correction factor might be an unpredictable function of model 
size and angle of attack. 


CONCLUDING REMARKS 


The solid-blockage interference for a doublet on the tunnel axis 
and the boundary interference for lifting wings in circular, rectangular, 
and two-dimensional slotted tunnels have been calculated by substituting 
an equivalent homogeneous boundary for the physical boundary of discrete 
slots. In the case of small wings, the results calculated with the 
assumption of homogeneity have been found to be consistent with those 
calculated for the discrete slots for as few as four slots in a circular 
tunnel. Furthermore, available experimental results for blockage inter- 
ference are consistent with the results of the present analysis. 

Through the concept of a wall restriction constant it is possible 
to reduce the interference due to all different slotted-wall configura- 
tions for the same tunnel shape to a single curve. Thus, the number of 
computations required to describe completely the interference for a 
given-shape slotted tunnel with different slot spacings and open ratios 
is greatly reduced. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., May 12, 1953- 
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Figure 1.- Schematic longitudinal and transverse cross-sectional views 
of a slotted tunnel with model installed, showing how air flows out 
and in through the slots as it passes the model. 
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Figure 2.- Sketch showing physical arrangement and nomenclature used in 
equation for potential difference across a thin, straight slotted wall. 



Figure 3* _ Sketch of the transformation of the section of the straight 
slotted wall between y = +l.d into an arc of a circular slotted wall . 



Open ratio, 
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Figure 4.- The function 
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Figure 5-- Comparison of the potential as given by equation ( 6 ) and by 
the approximate equation ( 7 ). d = n/b (8 slots); q = 7 . 



Figure 6.- Boundary-correction factors for lifting wings in circular 
tunnels with slotted walls. For thin slots, c = S. log e esc i rtr 
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slots c = log_ esc i jtr n , where N is the number of slots in 

rtNA e 2 0> 

horizontal wall and. r 0 is the open ratio of the horizontal walls. 
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slotted (horizontal) walls. 





